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ABSTRACT 

We conducted very long baseline interferometry (VLBI) observations of five radio-loud narrow-line 
Seyfert 1 (NLSl) galaxies in milliarcsecond resolutions at 1.7 GHz (A18 cm) using the Very Long 
Baseline Array (VLB A). Significant parsec-scale structures were revealed for three out of the five 
sources with high brightness temperature by direct imaging; this is convincing evidence for nonthermal 
jets. FBQS JI644-I-2619 with an inverted spectrum showed a prominent one-sided linear structure, 
indicating Doppler beaming with an intrinsic jet speed of > 0.74c. FBQS JI629-I-4007, also with 
an inverted spectrum, showed rapid flux variability, indicating Doppler beaming with an intrinsic 
jet speed of > 0.88c. Thus, we found convincing evidence that these two NLSls can generate at 
least mildly or highly relativistic jets, which may make them apparently radio loud even if they are 
intrinsically radio quiet. On the other hand, the other three NLSls had steep spectra and two of them 
showed significantly diffuse pc-scale structures, which were unlikely to be strongly beamed. Thus, 
some NLSls have the ability to generate jets strong enough to make them intrinsically radio loud 
without Doppler beaming. NLSls as a class show a number of extreme properties and radio-loud 
ones are very rare. We build on these radio results to understand that the central engines of radio- 
loud NLSls are essentially the same as that of other radio-loud AGNs in terms of the formation of 
nonthermal jets. 

Subject headings: galaxies: active — galaxies: jets — galaxies: Seyfert — radio continuum: galaxies 
— techniques: interferometric 



1. INTRODUCTION 

Active galactic nuclei (AGNs) are believed to be 
powered by gravitational energy released from accret- 
ing matter onto supermassive black holes (SMBHs) at 
the center of galaxies. Narrow-line Seyfert 1 galax- 
ies (NLSls) are a subclass of AGNs and are defined 
by the following optical spectral properties: (1) per- 
mitted lines are slightly broader than forbidden lines; 
(2) the fuU width at half-maximum (FWHM) of the 
II/3 emission line is <2000 km s^^; and (3) the flux 
ratio [O III]/H/3 <3, but exceptions are allowed when 
strong high-ionization iron lines such as [Fe VII] A6087 
and [Fe X] A6375, which are seldom seen in type- 2 
AGNs, are present (jOsterbrock fc Poggd 119851 : iPoggd 
I2OOOD . Other extreme properties of NLSls are strong 
permitted F e n emission lines in op tical and ultravio- 
let spectra (jBoroson fc Greeiil I1992D. rapid and large- 
amplitude X-ray variability (iLeighly '1999at ), a steep 
soft X-ray spectr um (Boiler et al. 1996; B randt et al.l 
ll997HGrupe et al.iri 998. 200 |), and a frequently observed 
blueshifted line profile dZamanov et al.l l2002t iBorosonI 
I2005t iLeighlv fc Moord I2004D . In terms of radio lumi- 
nosity (10^^-10^^ W Hz~^) and the compact source size 
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(< 300 pc; lUlvestad et a ll 119951: iMoranI |2000[ ). the ra- 
dio continuum properties of NLSls and normal Seyfert 1 
galaxies are not remarkably different. NLSls are gener- 
ally radio-quiet obje cts^. The fraction of radio-loud ob- 
jects i n NLSls is low (Stepan ian et al.l2003l : [Greene et al.l 
l2006HZhou eTal. 2006) ,' that is, ~ 7% (R > 10) and 
~ 2.5% (R > 100) IKomossa ct al. 2006a, see also 
IZhou fcWa^ 120021 : [Yuan et al...2008.) . which is signif- 
icantly lower than normal Seyf ert galaxies and quasars 
(- 10%-15%: [!^ezic et al.l[2002[ ). 

There is increasing evidence that the unusual 
optical/UV/X-ray properties of NLSls are physically 
related to the hi gh mass accretion rate s close to the 
Eddin gton limit (iBoroson fc Greed 119921: [Pounds et all 
119951 ISulentic et all 120001 : iBorosonI l2002t ) of the rela- 
tively low-mass black holes (~ 10^-10^ Mq), which are 
located at the l ower-mass end of the known popula- 
tion o f SMBHs ([Peterson et al.l [20001 : [Grupe fc Mathuil 
[2004 [Komossa fc Xu[ [20071 and references therein). 
Hence, NLSls are expected to be promising samples 
for understanding the nature of the extreme accre- 
tion onto rapidly growing SMBHs. The radio emis- 
sion of AGNs and X-ray binaries is suggested to 
be also r elated to the black hole mass and accre- 
tion rate (i Merloni ejJ] [2003|; [Heinz fc Sunvaev[ [2003[ : 
[Falcke etal . 2004). Many studies using large samples 
show a corr elation between radio loudness and the black 
hole i nass dLaoH [200d [Lacv et al.[ [2001[: [Dunlop et al. 



[20031: [McLure fc Jarvid [2004 [Metcalf fc Magliocchetti 



^ Radio loudness, _R, is conventionally defined as the ratio of the 
5-GHz radio to B-band flux densities, with a thresho ld of fl = 10 
sepa rating the radi o-loud and radio-quiet obje cts ( Visnovskv et al.1 
11992 : iStocke et aLiI 19921: IKellermann et al.|[l994 ') . 
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I2006t IWoo fc Urr^ |2002[ ) and an anti-correlatio n be- 
tween radio loudness and the accretion rate (iHol [20021: 
Lacy et alj 120011: iMaccarone et all 120031: iGreene et aLl 



20061) . These correlations predict that AGNs with high 



accretion rates and low-mass black holes such as NLSls 
should be radio quiet. In fact, as mentioned above, most 
NLSls have been detected as radio-quiet objects. In the 
analogy of X-ray binaries, NLSl s are at the high/soft 
state or at the very high state ([Mineshige et al.l 120001 : 
IWang fc Netzer "2003') . where radio emission is quenched 
(jMaccarone et al. 20031 . 

Radio- loud NLSl s are rare, but t hey do exist 
(iSiebert et al.l [Tool iGrupe et all [2OO0I: iZhou fc Wand 



2002t IZhou et al.l 



Komossa et al.l l2006allbl : lYuan et al.l I2008D 



200!" 



IWhalen et al.l l2Q0a 
The ra- 



dio powers of the radio- loud NLSls exceed those of 
the most luminous starburst galaxies, suggesting that 
strong AGN jet emissions are expected to be responsible 
for making these NLSls radio loud. Although the 
black hole masses of the radio- loud NLSls {^ 10^ -^^0) 
may be slightly larger than those of the radio-quiet 
NLSls, the radio-loud NLSls are still at a lower mass 
regime in AGNs and at a y ery high accretion rate 
regime (iKomossa et al.l l2006al see also IWhalen et al.l 
I2OO6I : lYuan et al.l 12008*): although such properties tend 
toward radio-quietness. One possible explanation of 
the existence of radio-loud NLSls is that the Doppler 
beaming effect on nonthermal jets influences radio 
loudness as well as some other radio-loud AGNs. 
However, most radio-loud NLSls show steep soft-excess 
X-ray spectra, whi ch are the same as those of the 
radio -quiet NLSls (jKomossa et al.l l2006al : lYuan et all 
l2008f ) and in contrast to the flat X-ray spectra due 
to the inverse Co mpton process in r elativistic jets in 
radio- loud quasars (jReeves et al.lllQQTI) . The steep X-ray 
spectra of radio-loud NLSls may be synchrotron tails 
originating froETL the relativistic jets (|Zhou et al.l 120071 : 
lYuan et al1l2008D , as seen in high-e nergy-peaked BL Lac 
objects (iPadoyani fc GioniinilllQQSD . On the other hand, 
IKomossa et al.l ( 2006aD indicated that many radio- loud 
NLSls show steep radio spectra where strong beaming is 
generally not expected, although very radio- loud NLSls 
may systematically show flat radio spectra ([Yuan et al.l 
|2008| ). Several radio-loud NLSls show physically unre- 
alistic, high brightness temperatures, derived from radio 
flux variability. This indicates the Doppler beaming 
effect if not int erstellar scattering (jZhou et al.l l2003t 
I Yuan et al.l 120081 ). However, the Doppler beaming effect 
on the jets of NLSls is not yet well established. 

Very long baseline interferometry (VLBI) is one of 
the most powerful tools for revealing the radio prop- 
erties in parsec (pc) scales by direct imaging at mil- 
liarcsecond (mas) resolutions. Arcsecond-resolution ob- 
servations have resolved the structures of only a few 
NLSls b ecause they are generally quite co mpact radio 
sources (l Ulvestad et"aLlll995t lMoranll2000[ ). VLBI ob- 
servations on a large number of radio-quiet and radio- 
loud NLSls are crucial for understanding the nature of 
the jet phenomena in these central engines. VLBI im- 
ages of only several objects have b een reported; that is , 
for radio-quiet NLSls: NGC 4395 (fVWobel et al.ll200l . 
MRK 766, AKN 564 (iLal et al.l [200l. and NGC 5506 
(jMid delberg et al.l [20041) a nd for radio-lou d NLSls: 
SDSS J094857.31-H002225.4 (iDoi et al.ll20"06a[ hereafter 



[Paper IT). J16290-f- 4007. J16333-I-4718, and B3 1702-1-457 
(|Gu fc Ghenl[2010h . SDSS J094857.31-h002225.4, which 
is the most radio-loud narrow-line object with R « 
2000, is su ggested to be Dop pler beamed by radio flux 
variability ([Zhou et al.l 120031 ) . and it has been revealed 
that Doppler-boosted jets are required to explain the 
directly measured high brightness temperatures of its 
very compact radio emissions seen in the VLBI im- 
ages ([Paper J ). The recent discovery of 7-ray emission 
from SDSS J094857.3 1-F002225.4 by the Fermi Gamma- 
Ray Space Telescope ([Abdo et al.ll2009al) and the subse- 
quent VLBI observations have established the presence 
of a relativistic jet in this extrenie nar row-line quasar 
([Abdo et al.ll2009bl : [ciroletti et al.l[20n[) . 

We conducted VLBI observations for five radio-loud 
NLSls at 1.7 and 8.4 GHz using the Very Long Baseline 
Array (VLB A) and the Japanese VLBI Network (JVN), 
respectively. The res ults of the JVN observati ons, whic h 
have been published ([Doi et al.ll2007L hereafter [Paper III ). 
showed no significant radio structure in these five ob- 
jects. In the present paper, we report the observations 
at 1.7 GHz, where synchrotron radio jet structures are 
generally expected to be more prominent than at higher 
frequencies. Based on the VLBA observations at 5 GHz 
for three radio-loud NLSls that are three of our five sam- 
ples (f or one source also at 2.3 and 8.4 GHz). [Gu fc ChenI 
([2010( ) reported that compact structures predominate ra- 
dio emissions. 

The present paper is structured as follows. The sample 
selection is shown in Section [2] Section [3| describes our 
observations and the data reduction procedures. In Sec- 
tion |4l we present the results and image analyses, and in 
Section[5l we discuss the observed pc-scale jet properties 
and their implications. Finally, the summary is presented 
in Section ini Throughout this paper, a flat cosmology is 
assumed, with Hq — 71 k m s~ ^ Mpc~^, JIm = 0.27, and 
Qa = 0.73 (ISoergel et al.l[2f)0l . 

2. SAMPLE 



IZhou fc Wa^ ([2OO2D found nine radio-loud NLSls 
from 205 NLSls listed in "A Catalogue of Quas ars a nd 
Active Nuclei: 10th ed. (fVeron-Cettv fc Veronl[200l ." 
Regarding these nine objects, we searched the Very Large 
Array (VLA) archival data obtained at 4.9 and 8.4 GHz 
with the VLA in the A-array configuration, which could 
provide radio positions with sufficient accuracy {'^ 0".l) 
for processing in VLBI correlators. We retrieved the 
data from the US National Radio Astronomy Observa- 
tory (NRAO) and analyzed them. As a result, we ob- 
tained radio positions for the five objects (Table[T|), which 
are our samples for the VL BA obser vations (this paper) 
and the JVN observations (iPaper III ). Thus, our sample 
selection processes were VLBI-oriented. No attempt was 
made to ensure sample completeness; a strict study on 
the general properties of the radio-loud NLSls is out of 
the scope of this paper. 

The properties of optical emission lines that are used 
to classify the five objects as NLSls in literatures are 
described in the Appendix; they are the FWHM of H/3 
of the broad component, the flux ratio [O III]/H/3, and 
^4570 J where -R4570 is the flux ratio of the Fe II mult iplets 
in the range 4434-4684A ([Veron-Cettv et al.|[200l[ ). 

3. OBSERVATIONS AND DATA REDUCTION 



Parsec-scale jets in radio-loud NLSls 



TABLE 1 

Radio-loud NLSl Samples for Our Observations 



Name 



(1) 



(2) 



^20cm -F'Scm log ij, 

(mjy) (mjy) 
(3) (4) (5) 



VLBl 



8.4GHz ^''SLj, 

(mJy) (erg s ^) 

(6) (7) 



RX J0806.6-I-7248 0.099 
FBQS J1629-I-4007 0.272 
RX J1633. 3-1-4718 0.116 
FBQS J1644-I-2619 0.145 
B3 1702+457 0.061 



50=^ 


26 


2.07 


5.0 


43.9 


12 


20 


2.15 


23.3 


45.3 


65 


35 


1.87 


22.1 


44.3 


91 


99 


2.97"^ 


150.3 


44.1 


119 


26 


1.02 


18.5 


44.0 



Note. — Col. (1) source name; Col. (2) redshift; Col. (3) 20- 
cm to tal radio flux density from the FIRST catalog llWhite et al.l 
119971): Col. (4) 6-cm tot al radio flux density from the GB6 cata- 
log I Cregorv et al.lll99B ): Col. (5) radio loudness, which was de- 
rived from the fc-corrected 6-cm ra dio flux density and t ^-band 
magnitude using the values listed in lVeron-Cettv &: Veroiil 1)20061 ) 
assuming an optical spectral index of —0.5; Col. (6) VLBI flu x 
density measured with the JVN at 8.4 GHz (A3. 6 cm; IPaper H) : 
Col. (7) ROSATX-ray lu minosi ty at 0.1-2.4 keV l IBade et al.lll995l : 
ILaurent-Muchlcis en et al.lll998h . 

"Flux density at 1.4 GHz from the NRAO VLA Sky Survey (NVSS; 
[Condon ct al. 199^ 

'^A much larger value compared to that in IZhou fc Wane! 
20021. mainly be cause they used my = 17.10 mag from 
Veron-Cettv &: Ve ron 200 1!, while we used my = 18.44 mag, 



on the basis of SPSS 4DR llAdelman-McCarthv et al.ii2006i) from 
I Veron-Cettv fc VeronI |[200^'l . 

The observations were performed using 10 antennas 
of the VLBA on 2006 February 07 and 10, while the 
JVN observations w ere condu cted on 2006 March 17 
and 26 and May 20 (jPaper III ): the difference between 
the dates of the VLBA and JVN observations was ~5- 
15 weeks. A left circular polarization was obtained at a 
center frequency of 1.667 GHz with a total bandwidth 
of 32 MHz. Because the targets were weak, we used a 
phase-referencing technique that involved fast switching 
of an antenna's pointing direction between a target and a 
reference calibrator in order to remove atmospheric fluc- 
tuation for increasing the coherent integration time. The 
switching-cycle period was usually 5 minutes or ^3 min- 
utes at low elevations^. 

Data reduction and imaging procedures were per- 
formed using the Astr o nomic al Image Processing Sys- 
tem (AIPS; [ Greisenl l2003l ) and Difmap software 
(|Shepherdlll997T ) , respectively. Amplitude calibration us- 
ing a priori gain values together with the system noise 
temperatures measured during the observations was ap- 
plied, and the calibration accuracy was typically 5%. In- 
appropriate earth-orientation parameters in the correla- 
tor were corrected. The geometrical correction of the 
antenna parallactic angles for the rotating feeds on the 
altitude-azimuth antennas in the position angle (P-A.) 
with respect to the source during the course of the ob- 
servation was performed. Ionospheric dispersive delay 
was corrected with a global ionospheric model using the 
task TECOR. Careful data inspection in time and spec- 
tral domains and data flagging were performed; there was 
no significant radio interference in the data. Bandpass 

® To make phase- refere ncing secure, more than one calibrator 
wa s used for each target llFomalont fc KopeikinI 120031 : IDoi et aLl 
1200 6b) during the JVN observation. However, this was not advan- 
tageous for the VLBA observation and significantly good results 
have been obtained without the secondary calibrators. 



correction was performed using the amplitude and phase 
data of a strong calibrator. The fringe-fitting solution 
for the reference calibrator was obtained and then ap- 
plied to the target ("phase- referencing"). We found that 
the phase-referenced visibilities of the targets had signif- 
icantly high signal-to-noise ratios for the self-calibration 
of the integrations for <30 minutes. We performed imag- 
ing with Difmap, which iteratively uses deconvolution 
and self-calibration algorithms for both amplitude and 
phase, starting from the phase-referenced data and a 
point source model. 

4. RESULTS 

We detected all five radio-loud NLSls. Images with 
significantly much higher dynamic ranges than those of 
the JVN were obtained. The resultant images are shown 
in Figure [T] the image parameters are listed in Table [2J 
At 1.7 GHz, all targets showed some structures in the 
VLBA images (see the Appendix for details on individ- 
ual objects), where as they w ere featureless in the JVN 
images at 8.4 GHz (jPaper IS ). The results of the mea- 
surements of peak intensities and flux densities are listed 
in Table 131 In this paper, we define peak intensity as the 
core in each image. 

4.1. Brightness Temperatures 

We calculated the brightness temperature for the core 
of each source. The brightness temperature in K of a 
component at the source frame is given by 



Tb = 1.8 X lO^(H-z) 



S. 



^ ^maj^Pmin 



(1) 



where z is the redshift, S^ is the observed flux density in 



mJy at frequency z/ in GHz, and 



^maj 



and (/)min in mas 



are the fitted FWHMs of the major and minor axes of the 
source size, respectively (jUlvestad et al.|[2"005[ ). Because 
the core in each source was unresolved, we adopted one- 
half of the beam sizes (i.e., 9/2; Tabled]) as the upper 
limits to the source sizes 0. The calculated brightness 
temperatures are listed in Table E] Tb > 10®'' K was 
derived for all sources. 

The geometrical mean of the synthesized beam sizes 
of the VLBA images was 4.6 x 10.2 mas and that of 
the image sensitivities was 0.144 mJy beam~^ (Table[2|), 
which was equivalent to the brightness temperature 2.0 x 
10^ K. 

4.2. Spectral Indices 

By comparing the flux densities (sum of the clean com- 
ponents), the spectral index a {Si, oc i''^") can be ob- 
tained (Table [3|) . The ranges of spatial frequency {uv- 
range), the ratios of the projected-baseline length to the 
wavelength, are ~ 1.0-50 MA at 1.7 GHz for the VLBA 
observations and ^ 4.0-65 MA for the JVN observa- 
tions at 8.4 GHz, which corresponds to the scales of 
4.1-210 mas and 3.2-52 mas, respectively. The similar 
uw-ranges and the proximity in time (difference of only 
~5-15 weeks) of the VLBA and JVN observations can 
minimize the systematic errors related to the spatial res- 
olution effects and flux variability in deriving a spectral 
index. Only RX J0806.6+7248'might have a possible 
structure larger than 52 mas, which may be resolved out 
in the JVN baselines at 8.4 GHz. In addition, the lower 
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sensitivity of the JVN images could derive the less clean 
components found above the image-noise levels of the 
images of diffuse emissions such as RX J0806. 6+7248. 
Hence, the derived spectral indices indicate at least a 
lower limit for RX J0806. 6+7248. 

To derive the spectral indices of the compact core com- 
ponents, we also express the spectral indices in terms of 
peak intensities, which defined as the strongest intensity, 
between 1.7 GHz and 8.4 GHz in the images. The resolu- 
tion effect is presumably negligible for these unresolved 
components. 

5. DISCUSSION 

We found significantly resolved pc-scale radio struc- 
tures for three NLSls and compact structures for two 
NLSls. The key results from our VLBA images are as 
follows: (1) detections of very high brightness temper- 
atures and linear structures of the pc-scale radio emis- 
sions, which provide evidence of nonthermal jets in the 
NLSls (Section 15. 1[) : (2) asymmetry of radio structures 
and core spectral indices that are very important infor- 
mation for evaluating the Doppler beaming effect on the 
jets (Section 15. 2|) : and (3) discovery of different radio 
properties of the five sources, which suggest why these 
NLSls are radio loud (Section l5.4|) . Detailed descriptions 
of individual objects are given in the Appendix. 

5.1. Strong Evidence for Nonthermal Jets 

The previous JVN observations at 8.4 GHz detected 
brightness temperatures of at least 10^'''^ K and high ra- 
dio luminosities, which cannot be attributed to any stel- 
lar process; radio emissions presumably originate from 
nonthermal jets. The VLBA observations have also re- 
vealed that radio components with very high brightness 
temperatures, which provide stronger constraints, i.e., 
Tb > 10^'^ K, were associated with the five NLSls. Sig- 
nificant pc-scale radio structures from the VLBA obser- 
vations were newly detected, providing stronger evidence 
for the presence of nonthermal jets. FBQS J1644+2619 
and B3 1702+457 clearly displayed linear structures. 
(RX J0806. 6+7248 also showed a linear structure in the 
east and west of its nuclear region.) 

RX J0806.6+7248, RX J1633.3+4718, and B3 
1702+457 are steep-spectrum sources {a < —0.5) with 
very high brightness temperatures. This is a strong ev- 
idence for a nonthermal process dominating the radio 
emissions. On the other hand, FBQS J1629+4007 and 
FBQS J1644+2619 showed inverted spectra (a > 0). 
The observed brightness temperatures were too high for 
free-free emissions, although a radio spectrum of not 
only nonthermal synchrotron emissions but also ther- 
mal free-free emissions could be inverted in an opti- 
cally thick regime^. We consider a putative source with 
having a diameter of lower than that of the VLBA 
beams radiating the free-free emissions at electron tem- 
peratures equal to the observed brightness tempera- 
tures and being optically thick at our observing fre- 
quency. Under these conditions, the X-ray luminosities 
at 0.1-2.4 keV of the bremsstrahlung emission would be 
^ 10^^ erg s^^, which cannot meet the observed X-ray 

''' For example, a thermal torus of at most T b ~ 10^ K such 
as the one inferred for type-2 Seyfert NGC 1068 IjGallimore et aO 
[T997l.[200D was imaged using the VLBA. 



luminosities ~ lO"*** erg s ^ (Table [It; these were de- 
tected in the ROSAT AW SkvSurvev ^Bade et all HM^ 
iLaurent-Muehleisen et al.lll998D and should include the 
contribution not only of bremsstrahlung emission but 
also the inverse Compton and synchrotron emissions. 
Hence, for all observed NLSls, we found strong evidence 
that the nonthermal process dominates the observed ra- 
dio emissions. 

Because the VLBI-correlated fluxes at 1.7 GHz were 
equivalent to much of the total fluxes in the VLA at 
1.4 GHz, we conclude that the radio emissions of all these 
radio-loud NLSls are dominated by the nonthermal syn- 
chrotron process in jets and that these pc-scale jets make 
the NLSls radio loud. 

5.2. Asymmetric Core-Jet Structures 

The AGN radio sources generally consist of steep- 
spectrum jets (a < 0) and a flat-/inverted-spectrum 
core (a > 0); the core is the location where the syn- 
chrotron emission becomes visible from an optically thick 
upstream to an optically thin downstream at an observ- 
ing frequency. The core flux contribution generally tends 
to a relative increase in the higher frequency radio im- 
ages; therefore, the JVN emissions should originate from 
cores. 

A structural asymmetry relative to the core was found 
in the limited dynamic ranges of our images, although 
only two components were seen in FBQS J1629+4007 
and RX J1633. 3+4718. The possible causes of the asym- 
metry are as follows: (1) the jets are intrinsically one- 
sided or more lossy (and hence radiative) on one side than 
the other; or (2) Doppler boosting effect on the relativis- 
tic jets; or (3) the jets are intrinsically symmetric but 
appear asymmetric owing to the stronger free-free ab- 
sorption toward the counter-jets due to an obscuring ion- 
ized gas intervening along our line of sight (|\yalker et al.l 
I2000ri Jones et all 120011 : iKameno et al.iboOlD . Both (2) 
and (3) require that the stronger side should be ap- 
proaching jets, while the weaker side should be reced- 
ing jets. In this paper, we discuss the possible pa- 
rameters of the jets in the NLSls on the assumption 
that the jet are intrinsically symmetric. The images 
of RX J0806.6+7248 and B3 1702+457 showed possible 
counter features, and thus indicating receding jets. How- 
ever, it should be noted that in some cases, the VLBI 
imaging processes probably generate a mimic counter 
co mponent near the c ore (e.g., imaging demonstrations 
by ICarilli et al.l Il994[ ) . We cannot be confident that 
the counter-jet features were genuine receding jets (the 
counter jets in RX J0806. 6+7248 may exist because a 
diffuse counter feature, which is misaligned with the lin- 
ear structure, was also seen at a small distance; the Ap- 
pendix. On the other hand, no counter component can 
be seen in FBQS J1644+2619, which has a significant 
linear jet structure in the limited image dynamic range; 
this result provides a reliable constraint for a receding jet. 
Hereafter, we consider FBQS J1644+2619 for discussing 
asymmetry. 

5.2.1. Doppler Beaming Effect Estimated from 
One-Sidedness 

In case Doppler boosting was responsible for the ob- 
served asymmetry of the jet structure, we can limit the 
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Fig. 1. — VLBA images of five radio-loud NLSls at 1.7 GHz. Source names are indicated at the upper left corner in each panel. All 
images were formed using natural weighting and were self-calibrated. Contour levels are separated by factors of \/2 beginning at three times 
the rms noise (Table [2ll. Negative and positive contours are shown as dashed and solid curves, respectively. Half-power beam sizes are given 
in the lower left corners. A large-scale map in panel (aO) was created by convolution using a restored beam of 15 mas (shown as an open 
circle) for easy recognition of the low surface brightness structures. The positions of the model-fitted components of FBQS J1629-(-4007 
and RX J1633. 3-1-4718 are indicated as points. 
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TABLE 2 
Parameters of VLBA Images 



Name 


(J 


t'niin ^ t'niaj 


RA. 


DR 




I 




(mjy beam^ 


■'■'■) (mas X mas) 


(°) 




(pc 


mas~^) 


(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


RX J0806.6+7248 


0.117 


4.7 X 8.4 


-20.5 


82 




1.8 


FBQS J1629+4007 


0.148 


4.1 X 10.7 


-10.1 


70 




4.1 


RX J1633. 3+4718 


0.114 


4.8 X 11.6 


35.3 


429 




2.1 


FBQS J1644+2619 


0.176 


4.5 X 9.2 


-4.1 


446 




2.5 


B3 1702+457 


0.166 


4.7 X 11.4 


-47.5 


357 




1.2 



Note. — Col. (1) source name; Col. (2) mis of image noise; Col. (3) FWHMs of minor and major axes of the synthesized beam; 
Col. (4) position angle of the beam major axis; Col. (5) image dynamic range, defined as the ratio of peak intensity (Table (3)1 to the rms 
of image noise; Col. (6) linear scale in pc corresponding to 1 mas at the distance to the source. 



TABLE 3 
Results of VLBA Observations at 1.7 GHz 



Source name 'j^^^ ■^i^VgHz log-Pl.TGHz logTg" 

(mJy beam-i) (mJy) (W Hz"!) (K) 



ai as 



M (3 (3) (4) (5) (6) (7) 

-0.51 ± 0.23 -0.94 ± 0.42 
+0.59 ± 0.19 +0.49 ± 0.26 
-0.73 ± 0.18 -0.57 ± 0.20 
+0.38 ± 0.16 +0.25 ± 0.16 
-0.84 ± 0.18 -0.90 ± 0.21 

Note. — Col. (1) source name; Col. (2) peak intensity; Col. (3) flux density of total CLEANed components; Col. (4) radio power at 
a rest frequency of 1.7 GHz, in which a fc-correction was applied assuming the spectral index shown in column (7); Col. (5) brightness 
temperature of the core component in the rest frame (Equation Q); Col. (6) spectral index derived from VLBA and JVN peak intensities; 
Col. (7) spectral index derived from VLBA and JVN flux densities. 



RX J0806.6+7248 


9.6 ± 0.5 


23.0 ± 1.2 


23.7 


>8.7 


FBQS J1629+4007 


10.4 ± 0.5 


10.5 ± 0.5 


24.2 


>8.9 


RX J1633. 3+4718 


48.9 ± 2.4 


55.5 ± 2.8 


24.3 


>9.3 


FBQS J1644+2619 


78.4 ± 3.9 


100.6 ± 5.0 


24.7 


>9.8 


B3 1702+457 


59.2 ± 3.0 


79.9 ± 4.0 


23.8 


>9.4 
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Parameters of 


TABLE 4 
THE DoppLER Beaming Effect 
FBQS J1644+2619 


FOR 


Source name 
(1) 


Comparison r k a /3min 
(2) (3) (4) (5) (6) 


*max 

(°) 

(7) 


FBQS J1644-I-2619 

t 


S^(core)/3fr >147 3 ai 0.74 
/^{jet)/3(T >37 2 as 0.77 


42 
39 



Note. — Col. (1) source name; Col. (2) measurements compared 
to derive the ratio of jet to counter jet. Si/ (core) means the core's 
flux density. 3cr means three times the rms image noise as an up- 
per limit for a counter jet. /^(jet) means the intensity of the jet 
emission at the position separated from the peak by the spatial res- 
olution along the jets; Col. (3) ratio of jet to counter jet; Col. (4) k 
value used for the calculation of the Doppler beaming effect, where 
q = k — a (Section l5.2.1ll : Col. (5) assumed spectral index (Tabled; 
Col. (6) lower limit of jet speed (Equation Q); Col. (7) upper limit 
of the angle of jet axis (Equation (|4]|). 

jet speed and inclination by comparing the approach- 
ing and receding jets on the assumption that the rest- 
frame radio spectra and intrinsic jet speeds of both the 
sides were the same. A bulk Lorentz factor is defined as 
r = (1 — /3^)~°-^, where /3 = v/c {v and c are the speeds 
of the source and light, respectively). The Doppler factor 
is defined as 

^-r(l-/3cos<i>)' ^^^ 

where $ is the jet viewing angle, which is the angle be- 
tween tlw_dJTection_ofthesource velocity and our line of 
sight (jGhisellini et al.lfl993l Appendix B). The frequency 
and the flux density at the rest frame, frost and 5'J^rcst: 
are observed as those at the observer frame, I'obs = (^i^rcst 
and SJ^obs = ^^SJ^rcst: respectively. Assuming the spec- 
tral index a, 5'°^§3 = S^-^'S'j;^^^ (iRvbicki fc LightmanI 
Il979| ). Hence, the fiux-density ratio of the approaching 
and receding jets is 



^rob.(/3.*- 



1 -I- /3 cos $ 
1 - ^ cos $ 



(3) 



where q — S — a; when considering the brightness ratio 
instead of the flu x-density ratio in coii tinuous jets q — 
2 — a is applied (|Ghisellini et al.|[i993[) . Therefore, we 
obtained 

^1/9 _ 1 

Using this inequality, we can estimate the lower limit of 
13 and the upper limit of $. 

Using the inequality, we derived the parameters of the 
Doppler effect for FBQS J1644-h2619 (Table g]) without 
considering the free-free absorption. We considered that 
the core was the first visible component in the approach- 
ing jet stream and that the fiux density of the counter 
component in the receding stream was less than three 
times the rms image noise (Table [2]). We found that 
/3 > 0.74, which is either a mildly or highly relativistic 
jet speed and that $ < 42°, which is a mild constraint for 
the inclination in this NLSl (Equation (|4])). For the case 
of the brightness ratio in continuous jets (q — 2 — a), 
we found a similar result of /? > 0.77 and $ < 39°. 



This result indicates that the central engine of an NLSl 
can generate nonthermal jets with quite a high speed. 
The radio loudness of this NLSl is potentially due to the 
Doppler beaming effect. 

We cannot obtain 5 by only considering the appar- 
ent jet structure. VLBI image monitoring will provide 
the apparent jet speed in order to solve both /3 and $ 
for obtaining d. Another method to obtain the con- 
straint 6 by imaging is to confirm the compactness of 
the source at a higher spatial resolution. For example, 
by VLBI observations, the very radio-loud {R « 2000) 
NLSl SDSS J094857.31-h002225.4 has highly relativistic 
jets viewed pole-on because of its compactness, which is 
indicative of an extremely high brightnes s temper ature 
requiring Doppler boosting of S > 2. 7-5. 5 (jPaper J ). 

5.2.2. Free-Free Absorption 

The possibility of the apparent jet asymmetry because 
of the free-free absorption due to the intervening ionized 
thermal plasma is discussed for FBQS J1644-(-2619. The 
optical depth of free-free absorption is 



Tff 



.235 X 10-^r: 



-1.35, 



c4L 



(5) 



where To (K) and rio (cm~'^) are the electron temper- 
ature and density of the absorbing gas, respectively, 
1/ is the frequency (GHz), and / is the path length 
in pc (M ezger fc HendersonI [19671 ) . The attenuation of 
more than a factor of 37 (Table 2]) that corresponds 
to Tff > 3.6 is required to make a putative counter-jet 
opaque. As a rough estimate, considering a path length 
equivalent to the length of the putative counter-jet of 
20 mas (50 pc), the fully ionized gas with To > 8000 K, 
and ly = (1 + ^:) x 1.7 GHz, a large amount of electrons of 
the column density of electrons Ne{= nj) > 10^^ cm~^ 
would be required in the nuclear region. Chandra X- 
ray observations of FBQS J1644-I-2619 showed that the 
fitted absorption column density iVu (natural gas) is al- 
most identical to the Galactic value 5.12 x 10^" cm~^ 
(jYuan et al.ll2008[ ). suggesting a very poor intervening 
medium toward the nucleus. Thus, such a large free-free 
absorber is unlikely in the FBQS J1644-I-2619 nucleus 
unless an extremely ionized fraction is assumed. Fur- 
thermore, the VLA image of FBQS J1644-I-2619 (Fig- 
ure [Ijd)) also seems to show an asymmetric structure, 
which is presumably related to the AGN activity (the 
Appendix). The large-scale asymmetry up to 30 kpc 
implies that the free-free absorption cannot be respon- 
sible for the apparent asymmetry, because such optically 
thick plasma would provide thermal radio emission of 
~ 4 Jy (Equation ([T])), which is much larger than the 
observed total flux density. Hence, we conclude that 
Doppler beaming should be preferentially efficient and 
it contributes to making FBQS J1644-I-2619 one-sided, 
unless the jets are intrinsically one-sided. The intriguing 
large structure of FBQS J1644+2619 is discussed later 
in Section [Ol 

5.3. Doppler Beaming Inferred from Flux Variability 

The Doppler factor can also be restricted by bright- 
ness temperature, which is derived from flux variabil- 
ity. For our sample, we checked the flux variation be- 
tween the NRAO VLA Sky Survey (NVSS: [Condon et al.l 
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TABLE 5 
Parameters of the Variability Doppler-beaming Effect for FBQS J 1629+4007 



Name 



(1) 



TFIHHT" 



-'1.4 
(mjy beam^^) 

(2) 



Obs. Date 



(3) 



•^1.4 
(mJy) 

(4) 



Obs. Date Variability log Tb 
(K) 
(5) (6) (7) 



(K) 
M W 



mm ^max Pmin 

(°) 

(10) (11) 



(W Hz-1) 

(12) (13) 



FBQS J1629+4007 

t 



11.97 ± 0.17 1994 Aug 21 9.0 ± 0.5 1995 Apr 16 5.6o- 13.0 11.2 4.0 14 0.88 <22.7 <0.74 

t t t t t t 12 2.1 28 0.63 <23.4 <1.43 



Note. — Col. (1) source name; Col.(2) 1.4-GHz peak flux density from FIRST; Col. (3) observation date of the FIRST data; Col. (4) 1.4- 
GHz peak flux density from NVSS; Col. (5) observation date of the NVSS data; Col. (6) significance of the detection of varia bility ; 
Col. (7) brightness temperature derived from variability; Col. (8) assumption of maximum intrinsic brightness temperature (Section 15.31 1; 
Col. (9) lower limit of Doppler factor; Col. (10) upper limit of viewing angle of jet axis; Col. (11) lower limit of jet speed; Col. (12) intrinsic 
radio power; Col. (13) intrinsic radio loudness. 

not only a nearly pole-on viewed orientation but also a 
relativistic bulk jet speed. The nuclei of these NLSls 
are expected to have the ability to accelerate the jets 
to the relativistic speed. In the projected distance, the 
asymmetry of FBQS J1644-I-2619 continues up to 30 kpc; 
kpc-scale asymmetry has also been seen in many blazars 
in th eir VLA images (jMurphv et al.l 119931: iKharb et all 
l201C)t) . The jet length 30 kpc corresponds to ~ 4 x 
10^° the Schwarzschild radius for FBQS J1644+2619jan 
estimated black hole mass of 8.4 x 10^ Mq; I Yuan et al.l 
2008). The l argest radio galaxi es have radio structures 
of > 1 Mpc (jPalmaet al.l 120001: iMachalski et~all [20081 . 
which corresponds to > 1 x 10^° the Schwarzschild ra- 
dius of a black hole mass of 10^ M©. The significance 
of FBQS J1644-h2619 is that at least some NLSl nuclei 
have the ability to be huge radio galaxies in black-hole- 
mass-normalized scales. A further study of the kpc-scale 
jets of FBQS J1644-f 2619 wiU be reported elsewhere. 

RX J0806.6-h7248 and B3 1702-^457, which are steep- 
spectrum sources with non-negligible diffuse structures, 
are unlikely to be strongly beamed; this indicates in- 
trinsic radio-loud sources without Doppler boosting. 
RX J0806.6-I-7248 shows jet bending at -20 pc from 
the core; this jet bending extends to ^^100 pc (Fig- 
ure [IJ . Such a structure resembles that of TeV gamma- 
ray blazar Mrk 501, which has inner jets strongly bend- 
ing at ~ 30 mas fr om the core an d conti nuing onto very 
diffuse outer jets. iGiroletti et al.l ()2004D suggested that 
highly relativistic bulk speeds and a progressive change 
in the jet direction (4-25°) with respect to the line of 
sight are possible and in agreement with the observed 
one-sided jet properties in the pc scale of Mrk 501. A 
quite small viewing angle is required to make such an 
intrinsically small bending ap pear large. The core spec- 
trum inverted up to 8.4 GHz ("G iroletti et al.ll200l) can 
attribute to a strong beaming effect on the spectral peak 
frequency I'obs = '^t'rest (Section l5.2.ip for pole-on viewed 
innermost jets. On the other hand, the inner jets of 
RX J0806. 6-1-7248 appear two-sided with respect to the 
core with a steep spectrum, and the extended jets origi- 
nate from the both sides of the two-sided structure (the 
Appendix); the extended jets of Mrk 501 distinctly orig- 
inate from the one-sized inner jets. It is difficult to at- 
tribute these radio properties of RX J0806. 6+7248 to 
beamed jets. 

Thus, the nuclei of some NLSls can generate jets that 
are strong enough to make them intrinsically radio loud. 
Consequently, either case (strongly or weakly beamed) 
suggests significant jet productivity in the radio-loud 



I1998D and the Faint I mages of the Radio Sky at Twenty 
centimeters (FIRST; iBecker et all I1995D at 1.4 GHz. 
Because of the resolution effect due to the smaller 
beam size of FIRST, we could find an intrinsic vari- 
ability only if the flux density in FIRST was larger 
than that in NVSS. Data for all sources except for 
RX J0806.6-I-7248 were available from both NVSS and 
FIRST. Only FBQS J1629-H4007 showed a FIRST flux 
density la rger than that in NVSS. Using the same tech- 
nique as iGhosh fc Punslvl (|2007D . we detected a flux 
variation in FBQS J1629-I-4007 with a significance of 
5.6cr; the brightness temperatures were derived to be 
j^qIS.o j^ (Table [5]). Based on the near equiparti- 
tion of ener gies between rad iating particles and a mag- 
netic field (jReadheadl 119941) or the catastrophic elec- 
tron energy loss due to an inver se Compton process 
([Kellermann fc Paulinv-Tothlll969l ). the brightness tem- 
perature of nonthermal radio emission is expected to be 
within a physically realistic upper limit of 10^^'^ K or 
10^^ K in the rest frame; the Doppler factor is expected 
to be 5 > 4.0 or 2.1 and the viewing angles should be 
14° or 28°, respectively. Thus, it is evident that the 
emitting source is a highly relativistic jet viewed pole- 
on. This source is presumably Doppler boosted and may 
be intrinsically in the radio-quiet population (Column 13 
of Table [S]). The variability Doppler factors for sev- 
eral objects have also been p reviously reporte d in a very 
radio-loud NLSl population ()Yuan et al.ll2008D . Our lim- 
ited dynamic-range VLBA image of FBQS J1629-h4007 
showed only a structure of two discrete components (the 
Appendix), which probably indicated a one-sided jet. It 
is important to conduct VLBI imaging more thoroughly 
and at more observing frequencies in the future. 

5.4. Nature of Radio-loud NLSl Radio Sources 

NLSls are expected to be radio-quiet AGNs on 
the basis of the possible dependence of radio loud- 
ness on both the black hol e mass and the accretion 
rate ()Boroson fc GreenI [l992l ). Even if all NLSls were 
intrinsically radio quiet, the Doppler beaming effect 
could explain why some NLSls appeared radio loud. 
Our VLBA study for the five radio-loud NLSls has 
shown both the cases that are likely and unlikely due 
to the Doppler beaming effect. Although they may 
be intrinsically radio quiet, FBQS J1629-I-4007 and 
FBQS J1644+2619 (inverted-spectrum sources) have 
shown evidence of a significant Doppler beaming effect, 
which presumably contributes to their radio loudness. 
We emphasize that the strong Doppler beaming requires 
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NLSls in our sample. NLSls as a class show a num- 
ber of extreme properties and radio-loud ones are very 
rare. However, we build on these radio results to under- 
stand that the central engines of radio-loud NLSls are 
essentially the same as that of other radio-loud AGNs in 
terms of the formation of nonthermal jets. 

NLSls show two types of X-ray continuum spectra. 
Normal (radio-quiet) NLSls show a prominent soft ex- 
cess plus a relativ e ly steep (F ~ 2 ) power-law hard tail 
(jHaba et al.l 120081: iLeighlyl '1999b). Both the compo- 
nents are proposed to originate from the accretion disks; 
the soft e xcess is thermal emission from the extremely 
slim disk (t Abramo wjcz et al.lll988HMineshige et al.ll2000l : 
I Wang fc N ctzcr 20Q3) and the power-law hard compo- 
nent originates from the scattered-up disk photons gen- 
erated bY_th£_Jiiverse_Coinpton process. On the other 
hand, lYuan et all (|2008[ ) recently indicated that while 
some radio- loud NLSls show X-ray spectra similar to 
normal NLSls, many radio-loud NLSls resemble the 
high-energy-peaked f iat-sp ectrum radio quasars (HFS- 
RQs; P erlman et aT| |1998[ ) in terms of the broadband 
spectral energy distribution (SED), which shows syn- 
chrotron spectra (due to beamed jets) peaked at UV/X- 
ray regimes, resulting in very steep X-ray spectra. 
FBQS J1629-H4007 and FBQS J1644-H2619, both of 
which are presumably Doppler b eamed sources i n our 
study, show HF SRQ-like SEDs (|Padovani et al.l [20021: 
lYuanet al.ll2008[ ). In particular, FBQS J1629-F4007 is 
the first confirmed HFSRQ for its modeled synchrotron 
peak at 2 X lO^fHz and its very steep X-ray spectrum. 
I Zhou et all ()2007t ) also suggested that the radio-loud 
NLSl J0324-h3410 was an NLSl-blazar composite with 
an HFSRQ-like SED. Although a small number of sources 
in our study are discussed here, the radio-loud NLSls 
with beamed jets may have jet-dominated HFSRQ-like 
SEDs. 

On the other hand, RX J1633. 3-^4718, which may be 
unlikely to be strongly Doppler beamed because of its 
steep spectrum, shows an X-ray spectrum with a soft ex- 
cess plu s a power-law hard t ail similar to several normal 
NLSls (jYuanet al.l [2001) . [Yuan et al.l (2010) showed 
this ultra-steep soft X-ray emission can be well fitted 
by an optically thick accretion disk, whose inferred pa- 
rameters (a black hole mass of ^ 3 x 10^ Mq and ac- 
cretion rate) are in good agreement with the indepen- 
dent e stimates ba. s ed o n an optical emission line spec- 
trum. lYuan et al.l (2008) discussed that the hard tail is 
flat (r = 1.37 ± 0.49) and may be an inverse Compton 
component originating fro m the jets, s i milar t o the flat- 
spectrum radio quasars. lAbdo et"al1 (|2009a[ ) also dis- 
cussed that the simultaneous optical/X-ray/7-ray SED 
of the very radio-loud NLSl SDSS J094857.31-I-002225.4 
showed a blue continuum attributed to the accretion disk 
and a hard X-ray spectrum attributed to the jet. Thus, 
the contribution of the disk could not be negligible even 
in the radio-loud NLSls. Perhaps, even the HFSRQ-like 
SED may be due to only the dominance of the thermal 
emission from the extremely slim disk. We would have 
to compare many SEDs and beamed-jet profiles by con- 
ducting more VLBI imaging observations for radio-loud 
NLSls in order to understand the origin of the X-ray 
spectra. 

Similar to NLSls, bro ad absorption line (BAL) quasars 
(jWevmann et al.lll991| ) are also considered to be highly 



accreting black h ole (with large r masses compared to 
NLSls ) systems (iBorosonI l2002f l and generally radio 
quiet (jBecker et al.l 120011 ) . The rarity of the BAL 
quasars as radio galaxies ([Gregg et al.l l2006l) is an im- 
portant similarity to NLSls . A VLBI study of BAL 
quasars (iJiang fc Wand [20031: [Montenegro- Monteset al.l 



[20091: iDoi et all [20091: iKunert-Bairaszewska et al.ll2010D 
is crucial for revealing their accreting and outflowing phe- 
nomena and determining their speed and orientation of 
the nonthermal jets by observing them. An understand- 
ing of the natures of the highly accreting systems is re- 
quired to explain both the subclasses of NLSls and BAL 
quasars. 

6. SUMMARY 

The VLB A observations at 1.7 GHz by direct imag- 
ing of pc-scale radio structure have provided convincing 
evidence of nonthermal jets in radio-loud NLSls. We 
resolved significant pc-scale structures of three NLSls; 
the pc-scale radio emissions provided most of the total 
fluxes found in the VLA resolutions in all the targets. 
The previously reported JVN images at 8.4 GHz enabled 
us to determine that our sample contained two inverted- 
and three steep-spectrum sources. FBQS J1644-I-2619, 
an inverted spectrum source, showed a remarkable one- 
sided pc-scale jet structure that indicated either a mildly 
or highly relativistic jet speed of /3 > 0.74 and an in- 
clination of < 42° on the basis of the framework of 
the Doppler beaming effect. FBQS J1629-I-4007, an 
inverted-spectrum source, showed compact radio struc- 
tures in the pc-scale and a significant flux variability 
between the data from NVSS and FIRST. The phys- 
ically unrealistic brightness temperature inferred from 
the variability requires a pole-on viewed relativistic jet 
of a Doppler factor oi 5 > 4.0, implying a jet speed of 
13 > 0.88 and an inclination of < 14°. The Doppler 
beaming effect is possibly responsible for the radio 
loudness in FBQS J1644-1-2619 and FBQS J1629-h4007. 
RX J0806.6+7248 and B3 1702+457, which are steep- 
spectrum sources with non-negligible diffuse structures 
in the pc scale, were unlikely to be strongly beamed; this 
indicates that they are intrinsically radio loud even with- 
out the beaming effect. We found both types of radio- 
loud NLSls with relativistically beamed jets and intrinsi- 
cally strong jets, as well as seen in other radio-loud AGN 
classes. 
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and Communications Technology. 

APPENDIX 
INDIVIDUAL OBJECTS 

We describe below the optical spectroscopic classification in the literature, the radio properties observed by our 
observations in terms of the pc-scale structure, and the radio spectrum of each NLSl, particularly from the viewpoint 
of the Doppler beaming effect. 

RX J0806. 6+7248 

The optical emission line properties of RX J0806.6+7248 are FWHM(H/3) =1060 km s~\ [O III]/H/3 < 3, the 
unambiguously presence of strong Fe II lines; thus, it is classified as an NLSl (jXu et al.lll999t l. Our VLBA image of 
RX J0806. 6+7248 at 1.7 GHz showed the most remarkable radio structure in our sample. Low-brightness emissions 
with temperatures 10^-10^ K extended northeast up to 100 pc from the intensity peak (Figure[TJaO)). The inner region 
of this source (Figure [ija)) displayed an east-west linear structure that consisted of four discrete components with high 
brightness temperatures (> 10* K), which were identified by visibility-based model fitting. We cannot completely rule 
out the possibility of a continuous linear structure. Such structures in this steep-spectrum source (Table |31) are direct 
evidence that nonthermal jets are responsible for the observed radio emission and for making this NLSl radio loud (a 
substantial fraction of the total fiux was not resolved out in the VLBA image; Tables [T] and |3|) . A two-sided structure 
with respect to the core can be seen in the linear structure. We note that diffuse emissions appear to emanate toward 
the north from both ends of the inner linear structure, which supports a genuine two-sidedness of the inner linear 
structure. The steep spectrum core and the two-sidedness are unlikely to be caused by a strong Doppler beaming 
effect in this NLSl, i.e., the jets are sub-relativistic or significantly inclined. 

FBQS J1629+4007 

The optical emission line properties of FBQS J1629-h4007 are FWHM(H/3) = 1238 ± 20 km s'^ [O III]/H/3 = 
439 ± 5/1563 ± 17, and -R4570 = 0.48 ± 0.02, which wer e analyz ed on the basis of S loan Digital Sky Survey (SDSS) 
DR3; thus, it is classified as an NLSl ()Zhou et al.ll2006t see also lWhalen et al.l l2006^. Our VLBA image at 1.7 GHz 
showed an apparently compact pc-scale structure, but visibility-based model fitting indicated two discrete components 
(~ 10 mJy and ^ 1 mJy) rather than a single elliptical Gaussian (Figure. [TJb)). The weaker component was seen at 
6.2a in a residual image after the subtraction of the stronger component, and was separated at P. A. = 110 ± 16° by 
2.7 ±0.6 mas (11 pc), which is more than half of the synt hesized beam. A V LBA image at 5 GHz also shows a slightly 
resolved component separated by 1.5 mas at P. A. = 122° ()Gu fc Chenll2010r) . The consistent position angles may imply 
the jet direction of this NLSl. Because the total flux was detected in the VLBA image (Tables [1] and [3]) , the radio 
emission of this NLSl was dominated by these pc-scale components with the high brightness temperatures (> 10* K). 
Although the jet feature was insufficient to discuss one/two-sidedness, the combination of the inverted radio spectrum 
and the very compact structure was suggestive of the Doppler beaming effect. We also show further evidence of 
Doppler beaming with the Doppler factor 5 > 4.0 (jet speed > 0.88c) from fiux variability in the NLSl in Section [5731 

RX J1633.3+4718 

The optical emission line properties of RX J1633.3-h4718 are FWHM(H^) = 909 ± 43 km s^^, [O III]/H/3 = 
902 ± 21/918 lb 8, i?457o = 1.02, which were analy z ed on the basis o f SDSS DR5; thus, it is classified as an NLSl 
(| Yuan et al.l 120081 see also IWisotzki fc Badd [l997t IZhou et al.ll2006[ ). Our VLBA image of RX J1633.3-H4718 at 
1.7 GHz showed an apparently compact pc-scale structure, but indicated two discrete components (^ 38 mJy and 
~ 18 mJy) rather than a single elliptical Gaussian in the visibility- based model fitting (Figure [TJc)). The weaker 
component was seen at 1 la in a residual image after the subtraction of the stronger component and was separated at 
P. A. — 79° ± 2° by 2.9 ± 0.1 mas (6.1 pc), about half of the synthesized b eam. A VLB A image at 5 GHz also showed 
a second component separated by 3.7 mas at P. A. — 79° from the core ()Gu fc Chenll2010ri . The consistent position 
angles may imply the jet direction of this NLSl. Because most of the total flux was detected in the VLBA image 
(Tables [U and [3]) , the radio emission of this NLSl was dominated by the pc-scale structure with very high brightness 
temperatures (> 10^ K). Unlike FBQS J1629+4007, this NLSl showed a steep spectrum, despite the very compact 
radio structure. Such radio properties cannot generally be attributed to strong Doppler beaming but were associated 
with the compact radio lobes with a sub-relativistic speed similar to tha t seen in the compact steep-spectrum objects, 
which are probably an AGN class of young radio galaxies ()0'Dealll998D . 

FBQS J1644+2619 

The optical emission line properties of FBQS J1644-I-2619 are FWHM(H/3) = 1507 ± 42 km s"\ [O III]/H/3 = 
120 + 4/1108 ± 22 , and Ra^j^ =0.75, which were analyzed on the basis of SDSS DR5; thus, it is classified as an NLSl 
()Yuan et al.ll2008t Iwhalen et al.ll2006[) . Our VLBA image of FBQS J1644-b2619 at 1.7 GHz showed a continuous linear 
pc-scale structure with a length of ~ 20 mas {^ 50 pc; Figure [Ijd)); the discrete components cannot be identified, and 
the width in the direction perpendicular to the linear structure cannot be resolved. No diffuse emission component 
was found even in a fairly high dynamic range (DR = 429; Table [5]); on the other hand, RX J0806. 6-1-7248 showed 
significant diffuse emissions even with a much lower dynamic range (DR = 82; Table [2]). The clear linear structure with 
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high brightness temperatures (> 10^ K) was a direct evidence that nonthermal jets were responsible for the observed 
radio emission in this NLSl and for making it radio loud (most of the total flux was detected in the VLBA image; 
Tables [Hand [3]). The one-sidedness together with the inverted spectrum was highly suggestive of the Doppler beaming 
effect. We determined the jet speed > 0.74c (Section F5.2.ip . The one-sidedness seems to be maintained up to the kpc 
scale, shown in the FIRST archival image (Figure [ijd)); this image was also shown bv IWhalen et al.l (|2006[ ). Using 
the AIPS task JMFIT, we identified two radio components separated by 11.7 arcsec, corresponding to ~30 kpc. The 
brighter component was almost unresolved and coincident with the position of the VLBI emissions (a cross symbol 
in the FIRST image), while the secondary component was slightly resolved. There is no optical counterpart for the 
secondary component in the NASA/IPAC Extragalactic Database. The P.A. (71°) and the structural elongation (72°) 
of the secondary component are consistent with the direction of the pc-scale structure (61°. 0-63°. 6). Hence, the 
secondary component was likely to physically relate to the nucleus of FBQS J1644-I-2619. 

B3 1702+457 

The optical emission line properties of B3 1702-^457 are FWHM(H/3) = 1040 km s~\ [O IIll/H/3 ^ 2.48, and 
i?4570 = 1-86; thus, it is classified as an NLSl ()Leighlvlll99"9bt IVeron-Cettv et al.|[200lt IWisotzki fc Badel [l997V Our 
VLBA image of B3 1702-1-457 showed an elongated structure in the southwest direction with a length of ~ 30 mas {^ 
35 pc; Figur e [He)); the elongat ed structure seemed to curve from the west to the south. The VLBA images at 5 GHz 
reported by iGu fc ChenI ()201 0f) showed a second weak component separated toward the west from an intensity peak 
in a much higher spatial resolution, which was consistent with our 1.7 GHz image around the core of the curved jets. 
Unlike FBQS J1644+2619 with a linear structure, the radio structure consisted of quite diffuse components (~ 10^ K), 
except for the intensity peak (> 10^'^ K). Most of the total flux was detected in the VLBA image (Tables [T] and |3]) . 
The elongated structure with high brightness temperatures in this steep-spectrum source (Table |3]) is evidence that 
nonthermal jets are responsible for the observed radio emission and for making this NLSl radio loud. A weak counter 
feature was apparently seen with respect to the core, but this feature may be a mimic (Section 15. 2p . We were not 
confident of two-sidedness in this NLSl. The steep-spectrum core and the reasonably diffuse jet are unlikely to be 
caused by the strong Doppler beaming effect in this NLSl, i.e., the jets are sub-relativistic or significantly inclined. 
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